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We find an order-disorder transition from crystals to disordered crystals for static packings of 
frictionless spheres. While the geometric indicators are mostly blind to the transition, disordered 
crystals already exhibit properties apart from crystals. The transition approaches the close packing 
of hard spheres, giving rise to the singularity of the close packing point. We evidence that both 
the transition and properties of disordered crystals are jointly determined by the structural orders 
and density. Near the transition, the elastic moduli and coordination number of disordered crystals 
show particular pressure dependence distinct from both crystals and jammed solids. 
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Order and disorder constitute two fundamental themes 
in condensed matter physics and materials science. Per¬ 
fect crystals, the epitome of order, provide an important 
starting point for understanding properties of solids [Ij. 
In contrast, materials such as glasses and granular assem¬ 
blies are highly disordered, which exhibit a set of univer¬ 
sal properties distinct from their crystalline counterparts 
[^-y- Expectedly, a crystal can evolve away from the 
perfect crystalline order and develop into an amorphous 


state when disorder is introduced [8|, [9j. While intensive 
efforts have been invested into properties of amorphous 
solids, the characterization of the regime between crystals 
and amorphous solids, especially the crossover between 
the physics of crystals and that of disordered solids, has 
not been carefully tackled, leaving the boundary between 
the two extremes of order and disorder vague. 

Intuitively, one may distinguish disordered solids from 
crystals based on the structural orders and simply clas¬ 
sify solids with high structural orders to crystals. This 
has been proven infeasible because some solids with ex¬ 
tremely high structural orders could exhibit features of 
disordered solids [31 [llj : whereas those with consider¬ 
ably low structural orders may respond like crystals [111. 
Moreover, a recent experiment showed that a glass could 
exhibit low-temperature thermodynamic properties like 
polycrystals when being compressed to high pressures 
[l^ . This suggests that in addition to the structural or¬ 
ders there exist other possible controlling parameters of 
the manifestation of disordered solids, while the density 
is a candidate. However, how the density works together 
with the structural orders to determine properties of dis¬ 
ordered solids is still an open question. 

Bearing in mind these puzzles, we numerically inves¬ 
tigate the evolution from perfect crystals to disordered 
solids, by successively tuning the particle-size polydis- 
persity. In the similar framework, previous simulations 
have characterized the structural amorphisation at a suf¬ 
ficiently large polydispersity Here we focus on 


another unknown order-disorder transition at a rather 
small polydispersity from crystals to disordered crystals, 
namely solids with extremely high crystalline order in 
structure but mechanical and vibrational properties re¬ 
sembling disordered solids. While the structural orders 
i i i m are insensitive to this transition, multiple 
quantities undergo apparent changes. The critical poly¬ 
dispersity of the transition r]c is scaled linearly with the 
packing fraction distance from the close packing of hard 
spheres, (j) — ^cp- Therefore, the close packing behaves 
like a singular point where infinitesimally small polydis¬ 
persity turns the crystal into a disordered crystal M- 
The significance of this transition is further manifested 
by some unknown physics, which unveil important as¬ 
pects of solids and answers the questions raised above. 

We start with a perfect crystal, i.e., triangular lattice 
in two dimensions (2D) and face-centered cubic lattice 
in three dimensions (3D). Periodic boundary conditions 
are applied in all directions. All particles have the same 
mass m and interact via the potential 



where is the separation between particles i and j, cr^ 
is the sum of their radii, and 0(x) is the Heaviside func¬ 
tion. It is a valid model of experimental systems such 
as granular and colloidal solids [IM3. The disorder is 
continuously introduced by tuning the particle-size poly¬ 
dispersity. Particle i is assigned a random number Xi 
uniformly distributed in [—0.5,0.5], which leads to the 
particle diameter ai = {1 + Xiri)a with rj the polydis¬ 
persity. We increase rj from 0 successively by a small 
step size Arj G [1.6 x 10“®,8 x 10“^], with smaller Arj 
applied near the close packing. After each change of rj, 
we rescale the average particle radius to maintain a fixed 
packing fraction. The system is then relaxed to the local 
potential energy minimum by the fast inertial relaxation 
engine minimization method (3 |. To obtain the normal 
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FIG. 1: Phase diagram with two order-disorder transitions 
in the parametric space of the polydispersity rj and packing 
fraction The circles with a linear fit rjc ~ (solid 

line) mark the transition from crystals to disordered crystals. 
The transition labeled by the squares signals the structural 
amorphisation from disordered crystals to amorphous solids 
(see the Supplemental Material for the definition). The 
dashed line is to guide the eye. The color contour shows 
log(l — ^e) with Te the bond orientational order. The dot- 
dashed line marks = 0.95, to the left of which all states 
have a high crystalline order Tg > 0.95. 


modes of vibration, we diagonalize the dynamical ma¬ 
trix using ARPACK 2]|. The mass, energy, and length 
are in units of the particle mass m, characteristic energy 
scale e, and average particle diameter a. We show here 
results of 2D packings with harmonic repulsion (a = 2). 
If not specified, the number of particles is A^ = 1024. 
We have verified that our major findings are valid for 
packings with Hertzian repulsion (a = 5/2), in both 2D 
and 3D, and with other polydispersity distributions, e.g., 
Gaussian. 

Figure [1] shows a unified phase diagram of solids over 
the entire spectrum of disorder with both the transi¬ 
tions from crystals to disordered crystals and from dis¬ 
ordered crystals to amorphous solids. Here we denote 
disordered solids with strong structural amorphisation 
as amorphous solids to distinguish them from disordered 
crystals. 

When a crystal is driven progressively into a disor¬ 
dered crystal, the contact network distorts successively 
and is eventually destroyed by some local contact break¬ 
ing. Consequently, the average coordination number z, 
i.e., the average number of particles with which a given 
particle interacts, drops below 6, as shown in Fig. [2](a). 
The contact breaking happens randomly in space, re¬ 
sulting in the spatially heterogeneous disorder, which is 
one of the most important features of disordered solids 

We thus propose Sz = — z)^ 

as the order parameter to characterize the strength of 
disorder, where Zi is the coordination number of parti¬ 
cle i. As shown in Fig. HKb), 6z increases quickly from 
zero at some polydispersity, signaling the transition from 


crystals to disordered crystals. 

To unambiguously determine the transition point, we 
calculate the susceptibility y = A^[(((5z)^) — where 

(.) denotes the average over 1000 distinct realizations un¬ 
der the same macroscopic conditions. The susceptibility 
method is superior in locating phase transition points 
Hi [ 2 ^ [ 2 ^ [^ . As shown in Fig. He), there is a peak 
in whose location rjc is defined here as the transition 
point from crystals to disordered crystals. As illustrated 
by the circles in Fig.[l] rje ^ 4> — (t>cp, where i/cp = 
is the packing fraction of close-packed hard disks. The 
close packing point is thus singular, because infinitesi¬ 
mally small polydispersity will trigger the transition. 

Interestingly, when y is plotted against r]/{(p — /icp) as 
in Fig. Hd), all curves collapse nicely, suggesting that 
ri/{(j) — (j)cp) be a more meaningful parameter in control 
of the transition. 

Figures[2](e)-(h) show the rj evolution of typical proper¬ 
ties characterizing disordered solids, including the elastic 
moduli and nonafBnity upon deformation. All quanti¬ 
ties change remarkably across rje- In the crystal regime, 
both the bulk modulus B and shear modulus G remain 
mostly constant in rj. Meanwhile, both the compression 
and shear deformations are affine with the corresponding 
nonafSnity /4c « 0 and ~ 0 (see Ref. [H for the defini¬ 
tion of /ic and /is). When rj > rje, B and G decrease, while 
/ic and /is increase, all at once. These changes strongly 
verify the validity and robustness of the transition. 

Figures[2](e)-(h) also demonstrate that disordered crys¬ 
tals at higher packing fractions and larger polydispersi- 
ties (stronger structural disorder) can have quantitatively 
similar mechanical properties to those at lower packing 
fractions and smaller polydispersities (weaker structural 
disorder). In addition to the transition at t/c, properties 
of disordered crystals seem to be also jointly determined 
by the structural orders and packing fraction in the form 
of rj/{4> — tpep)- This provides us with some clues to un¬ 
derstand the puzzle why compressed glasses can behave 
like polycrystals 0. When a glass is compressed, the 
ratio of the structural disorder to the density is smaller 
and approaches the value of polycrystals, which pushes 
properties of the glass closer to polycrystals. Therefore, 
it is the interplay between the structural orders and the 
density that determines the performance of a solid. To 
claim either of them to be deterministic is partial. 

Figure |3] manifest further the importance of disorder to 
disordered crystals from vibrational properties. One of 
the most special vibrational features of disordered solids 
is the boson peak, i.e., the peak in the reduced density of 
vibrational states with w the frequency and 

d the dimension of space (1.1^ . Figure Ha) shows 

D{(jj)/oj for disordered crystals at </ = 0.91. To smooth 
out the plane wave-like peaks due to finite size effect, we 
avera ge D jijj) at different system sizes from N = 256 to 
1024 [33 • With increasing r/, the boson peak (the first 
peak at low frequencies) gradually rises and moves to 
















3 



FIG. 2: (a)-(c) Polydispersity evolutions of the average coordination number z, fluctuation of the coordination number 5z as 
the order parameter, and susceptibility of the order parameter x- (d) Scaling collapse of all curves in (c) when x is plotted 
against ri/{(j> — 4>cp)- (e)-(h) Polydispersity evolution of the bulk modulus B, shear modulus G, nonafflnity of the compression 
deformation fj,c, and nonafflnity of the shear deformation fis- 


lower frequencies, consistent with the argument that the 
boson peak is correlated with the structural disorder 
[ 3 ^ . There are two other peaks of van Hove singularities 
at higher frequencies, whose presence indicates that the 
solids are still pretty crystalline in structure and possess 
hybridized characters of crystals and disordered solids. 

By plotting ZI(wbp)/wbp against 77 in Fig. [3jb) with 
wbp the boson peak frequency, we estimate below what 
value of 77 the boson peak disappears. Owing to the 
very strong finite size effect close to r]c, we cannot obtain 
smooth enough D{uj) to resolve the boson peak. By ex¬ 
trapolating the roughly linear part of the low 77 data, we 
find that D{u!Bp)/tOBP hits the Debye level at rjBp ~ rjc- 
As shown in Fig. me), rjBp ~ 77c over a wide range of 
packing fractions. The formation of the boson peak is 
thus another evidence to distinguish disordered crystals 
from crystals. 


Recently, it was proposed that solids spanning the en¬ 
tire spectrum of disorder could be described by either 
the physics of jamming or the physics of crystals till. 
For jammed packings of frictionless spheres, both the ex¬ 
cess coordination number z — zjso and the ratio of the 
shear modulus to the bulk modulus G/B are scaled well 
with the pressure P 38-41|: z — Ziso ^ pi/2(“-i) and 
G/B ^ with Ziso = 2d the isostatic value; 

while for crystals both z and G/B are independent of 
the pressure. With a more accurate control of structural 
orders, we check how the two types of physics evolve to 
each other. 


persity. Figured shows the pressure dependence of G/B 
and z — Ziso at different 77. When 77 is large, the jamming 
scalings are recovered over the whole range of pressures 
studied. When 77 is small, there is a clear transition from 
the high pressure crystal scalings to the low pressure scal¬ 
ings disobeying both the physics of crystals and that of 
jamming: G/B ^ pO-73 z —ziso p 0 . 86 _ transi¬ 
tion is right at the transition from crystals to disordered 
crystals. At intermediate 77, the three types of scalings 
are all present, with the newly reported scalings sitting 
between those of crystals and jamming. Therefore, dis¬ 
ordered crystals close to the transition at 77 c comprise a 
third family of solids complying with the physics other 
than crystals and jamming, whose origin and underlying 
physics are interesting issues to explore. 

The finding and characterization of the transition from 
crystals to disordered crystals reveals some unknown fea¬ 
tures: (i) The close packing is singular in terms of the 
transition, implying that it is the only rigid packing of 
hard spheres satisfying the physics of crystals; (ii) the 
structural orders and density interplay to determine the 
transition and properties of disordered solids; (hi) disor¬ 
dered crystals near the transition exhibit unique pressure 
scalings apart from crystals and jamming. Here we man¬ 
ifest that our knowledge about solids, even about seem¬ 
ingly crystalline solids, is still rather incomplete. Follow¬ 
up studies, especially to determine the nature of the tran¬ 
sition, are necessary to have a deeper understanding of 
the phenomena reported here. 


We start from a packing aX </> = 0.92 and quasistatically 
decrease the packing fraction (pressure) at fixed polydis- 


Our major findings may not be limited to ather- 
mal systems with repulsions. For systems with long- 
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FIG. 3: (a) Reduced density of vibrational states D(uj)/uj 
at (j> = 0.91 and different polydispersities. From the left to 
the right, the three peaks are respectively the boson peak and 
two van Hove singularities, (b) Polydispersity evolution of the 
strength of the boson peak D{ljJbp)/ujbp from (a). The solid 
line is a linear fit to the low rj data. It hits the Debye level 
labeled by the horizontal dashed line at t^bp. (c) Correlation 
between tjbp and rjc- The data points are calculated at (f> = 
0.907, 0.908, 0.91, 0.92, and 0.94 in the ascendent order of rjc- 
The dashed line shows r^BP = r)c- 
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FIG. 4: Pressure evolutions of (a) the ratio of the shear mod¬ 
ulus to the bulk modulus G/B and (b) excess coordination 
number z — 2 iso. The horizontal dotted lines show the crys¬ 
tal behavior. The solid (dashed) lines are power law fits to 
the data: G/B ~ P° ® {G/B ~ and 2 - 2iso -- ® 


range attractions, we expect similar results, based on 
the observation that strong repulsive interactions gov¬ 
ern the behaviors of the packings while attractions act 
as perturbations [i^. Since the effective interaction of 
some widely studied soft colloidal particles (e.g., poly-N- 
isopropylacrylamide) can be well fitted to the repulsions 
used here [l3, [3 : observations are also highly rele¬ 

vant to colloidal experiments, but the temperature effects 
should be also included. Our preliminary results show 
that the transition from crystals to disordered crystals 
shifts continuously to lower values of the polydispersity 
with increasing temperature because the thermal 
fluctuation is another factor in addition to the polydisop- 
ersity to induce frustrations. Therefore, the athermal 
case discussed here is not singular from thermal systems 
and the physics can be generalized to thermal systems 
well below the melting temperature, which require fur¬ 
ther studies and verification of colloidal experiments. 
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SUPPLEMENTAL MATERIAL 


Section I. Structural amorphisation from disordered 
crystals to amorphous solids 


As shown in the main text, a crystal undergoes a hid¬ 
den order-disorder transition into the disordered crystal 
phase at a very small particle-size polydispersity. In this 
section, we focus on the order-disorder transition at a 
sufficiently large polydispersity labeled by the squares in 
Fig. 1 of the main text, namely the structural amorphi¬ 
sation [H, 0. 

Here we calculate two widely used geometric order pa- 
rameterSjthe bond orientational order and translational 
order [3,0] . The bond orientational order of particle j is 
given by 


= 


■E' 


, 62 ^-j 


( 2 ) 


where rij is the number of nearest neighbors of particle 
j, and 9jk is the angle of the bond between particle j 
and its neighbor k with respect to the x axis. The global 
bond orientational order is the average over all particles: 
4'6 = The translational order is defined as 




1 

N 


N 

i=i 


(3) 


where G is any first shell reciprocal lattice vector of a 
hexagonally close packed system. For a perfect hexag¬ 
onal lattice, both 'I'e and T't are equal to one, while 
they decrease to small values for highly disordered solids. 
In order to unambiguously determine the amorphisation 
transition point, we calculate the order parameter sus¬ 
ceptibilities 


Ae.T = Nd'i'l j,) - (4) 

where (.) denotes the average over 1000 realizations under 
the same macroscopic conditions. 

Figure 3S1 shows the polydispersity evolution of the 
pressure P, potential energy per particle U, order pa¬ 
rameters 'I'e and 'Ft, and susceptibilities xq and xt cal¬ 
culated at various packing fractions. With increasing 
the polydispersity ry, P and U reach their peak values 
at r]a ~ 0.4. Meanwhile, both fFe and dfy show a fast de¬ 
cay. The behaviors of ^Fg and fF-T across rja resemble that 
of the two-dimensional melting [J|. As demonstrated by 
Figs. gsTc) and 930, the transition points located by 
the peaks in Xeiv) and XTiv) are identical, which signal 
a one-step transition from disordered crystals to amor¬ 
phous solids. The transition of a (disordered) crystal into 
a glass, driven by quenched disorder, has recently been 
observed in experimental systems of colloidal suspensions 

i- 
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FIG. S5: Polydispersity dependence of (a) pressure P, (b) bond orientational order parameter (c) susceptibility of the 
bond orientational order X6, (d) potential energy per particle U, (e) translational order parameter 'i’r, and (f) susceptibility of 
the translational order xt- Both P and U show a peak when the system undergoes the structural failure. At the same point, 
^<6 and show a fast decay and their susceptibilities also exhibit a peak. The amorphisation transition signaled by the peak 
in X6 is shown by the squares in Fig. 1 of the main text. 


Note that the critical packing fraction rja of the amor- 
phisation transition is almost independent of the packing 
fraction, in contrast to the linear scaling of the transi¬ 
tion from crystals to disordered crystals with respect to 
the packing fraction reported in the main text. More¬ 
over, Figs. and show that the geometric order 

parameters as a function of the polydispersity are also 
almost independent of the packing fraction, which are 
important evidences supporting the generality of Fig. 1 
of the main text. If the polydispersity is replaced with 
either of the geometric order parameters or even other 
structural order parameters, Fig. 1 of the main text is 
still a valid phase diagram. 

Section II. Alternate route probing the transition 
from crystals to disordered crystals 

To demonstrate that the hidden order-disorder tran¬ 
sition reported in the main text is robust and indepen¬ 
dent of the route, we verify in this section that the same 
transition occurs on the route of decompression at fixed 
particle-size polydispersity. In Fig. 1 of the main text, 
this route is perpendicular to that of varying the polydis¬ 
persity at fixed packing fraction to obtain the transition 
line labeled by the circles. We start from crystal states 
at high packing fractions and decompress them at fixed 
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FIG. S6: Susceptibility of the order parameter x as a func¬ 
tion of the packing fraction (j) obtained by decompressing the 
conhgurations from 4> = 0.92 at fixed polydispersities rj. The 
transition points signaled by the peaks in x(0) h® perfectly 
on the transition line obtained on the route of varying rj at 
constant 4>- 


polydispersity by a small decrement of the packing frac¬ 
tion A(^ € [3.625 X 10“®,2 x 10“^], with smaller Af/> for 
smaller polydispersity. Potential energy minimization is 
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FIG. S7: Over a wide range of packing fractions, the de¬ 
viations from perfect crystalline order of bond orientational 
order parameter 1 — ^<6 (filled points) and transitional order 
parameters 1 — (open points) are plotted as functions of 
rj. The vertical dashed lines label the critical polydispersities 
of the transition. The solid lines have a slope of 2, indicating 
the power law scalings in the crystal regime. 


FIG. S8: The departure of the vibrational modes from those 
of the perfect crystal as a function of the polydispersity 
r}. The vertical dashed lines label the critical polydispersi¬ 
ties r]c of the transition from crystals to disordered crystals. 
Interestingly, the critical A„ at rjc are constant in packing 
fractions, as illustrated by the horizontal dashed line. 


performed after each decompression step. The packing 
fraction evolution of the order parameter susceptibility 
X defined in the main text is shown in Fig. ^ There is 
also a peak in x{4‘)j signaling the hidden order-disorder 
transition. We have verified (not shown) that the tran¬ 
sition points obtained here lie perfectly on the transition 
line (circles) shown in Fig. 1 of the main text. 


Section III. Signs from the geometric order 
parameters across the hidden order-disorder 
transition 


Seen from Figs. 93b) and 93®) ^ there is no observ¬ 
able sign of the transition from crystals to disordered 
crystals. However, we still expect to see some changes 
of the geometric order parameters across the transition. 
The idea is that the mechanical network breaks through 
the transition, which nontrivially results in the change of 
the elastic properties (e.g., the bulk modulus B and shear 
modulus G shown in Figs. 2(e)-(h) of the main text) and 
should also be reflected in the geometric response. We 
then plot 1 — 'Fe and 1 — IFt as a function of the poly¬ 
dispersity and focus on the vicinity of the transition in 
Fig. 911 In the regime of crystals, both 1 — and 1 — iFr 
are scaled well with rj^. Weak deviations from the power- 
law scalings can be observed after the crystals transit to 
disordered crystals. 


Section IV. Another signature of the transition from 
crystals to disordered crystals from the mode 
analysis 


Normal modes of vibration are the fundamentals to 
understanding properties of solids. In a finite-size crys¬ 
tal, the density of vibrational states D{uj) is composed 
of a set of (5-functions, with each 5-function contain¬ 
ing degenerate modes. In a weakly disordered crystal, 
the mode degeneracy is broken and the 5-functions are 
broadened with a finite width positively correlated with 
the strength of disorder. Therefore, only in the pres¬ 
ence of sufficiently strong disorder and/or for sufficiently 
large systems, D{uj) can be smoothed out. Constrained 
by the computational power, we are unable to study 
large enough systems to reliably resolve the boson peak 
at small polydispersities right above the transition from 
crystals to disordered crystals. Here we introduce an 
alternate characterization of the mode evolution with in¬ 
creasing polydispersity ry, which bypasses the issue of 
finite-size effect but shows an interesting and robust fea¬ 
ture of the transition from the perspective of vibrational 
modes. 

We directly trace the evolution of the mode frequen¬ 
cies in reference with those of a perfect crystal, which is 
calculated as 
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where d is the dimension of space, LOi and wf are the 
eigenfrequencies of the i-th mode of the solid with given 
polydispersity and of the perfect crystal, and (.) denotes 




















the ensemble average. Due to the periodic boundary con¬ 
ditions, there are d zero-frequency modes, so the total 
number of nontrivial modes is dN — d. All the frequencies 
are sorted in the ascending order. At sufficiently small 77 , 
the i-th mode of a slightly deformed crystal only slightly 
deviates from the f-th mode of the perfect crystal, so that 
is an exact evaluation of the mode deviation from a 
perfect crystal. At larger i] when mode interchanges take 
place, the i-th mode of the disordered solid may not be 
directly evolved from the z-th mode of the perfect crys¬ 
tal any more, so the expected one-to-one correspondence 
breaks. However, we still expect to be a good quanti¬ 
tative calculation of the deviation from perfect crystals. 
As shown in Fig. 33 A;^ increases with rj. In the crys¬ 
tal regime, A,^ ~ rj^. Near the transition from crystals 
to disordered crystals defined by the peak in the order 
parameter susceptibility, A,^ grows abruptly. Interest¬ 
ingly, A^j at the transition is independent of the packing 
fraction, as shown by the intersections between the data 
curves and vertical lines illustrated by the horizontal line 
in Fig. 33 This result is suggestive of a universal fea¬ 


ture of the transition from crystals to disordered crystals. 
In combination with the observation that the transition 
at rjc may also signal the emergence of the boson peak, 
Ri 1.6 X 10“® at rjc then sets a critical amount of devi¬ 
ation from a perfect crystal, above which the boson peak 
starts to appear. However, the origin of this particular 
value of Ai^ is unknown at the moment, which deserves 
further investigations. 
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